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ABSTRACT 

Context. Multi-wavelength (X-ray to radio) monitoring of Young Stellar Objects (YSOs) can provide important information about 
physical processes at the stellar surface, in the stellar corona, and/or in the inner circumstellar disk regions. While coronal processes 
should mainly cause variations in the X-ray and radio bands, accretion processes may be traced by time-correlated variability in the 
X-ray and optical/infrared bands. Several multi-wavelength studies have been successfully performed for field stars and ~ 1 - 10 Myr 
old T Tauri stars, but so far no such study succeeded in detecting simultaneous X-ray to radio variability in extremely young objects 
like class I and class protostars. 

Aims. Here we present the first simultaneous X-ray, radio, near-infrared, and optical monitoring of YSOs, targeting the Coronet clus- 
ter in the Corona Australis star-forming region, which harbors at least one class protostar, several class I objects, numerous T Tauri 
stars, and a few Herbig AeBe stars. 

Methods. In August 2005, we obtained five epochs of Chandra X-ray observations on nearly successive days accompanied by simul- 
taneous radio observations at the NRAO Very Large Array during four epochs, as well as by simultaneous optical and near-infrared 
observations from ground-based telescopes in Chile and South Africa. 

Results. Seven objects are detected simultaneously in the X-ray, radio, and optical/infrared bands; they constitute our core sample. 
While most of these sources exhibit clear variability in the X-ray regime and several also display optical/infrared variability, none 
of them shows significant radio variability on the timescales probed. We also do not find any case of clearly time-correlated opti- 
cal/infrared and X-ray variability. Remarkable intra-band variability is found for the class I protostar IRS 5 which shows much lower 
radio fluxes than in previous observations, and the Herbig Ae star R CrA, which displays enhanced X-ray emission during the last 
two epochs, but no time-correlated variations are seen for these objects in the other bands. The two components of S CrA vary nearly 
synchronously in the / band. 

Conclusions. The absence of time-correlated multi-wavelength variability suggests that there is no direct link between the X-ray and 
optical/infrared emission and supports the notion that accretion is not an important source for the X-ray emission of these YSOs. No 
significant radio variability was found on timescales of days. 

Key words, stars: pre-main sequence - stars: individual: R CrA, S CrA - radio continuum: stars - X-rays: stars - Infrared: stars 



1. Introduction 

Before low-mass young stars reach the main sequence they ma- 
ture throu gh a serie s of evo lutionary stages te rmed class to 
class III (Lad all 19871 and e.g. lAndre et aTll2000h . Named 'Young 
Stellar Objects' (YSOs) throughout this development, in the ear- 
liest stages they are referred to as 'protostars'. A class pro- 
tostar, observed as cold submillimeter source still has to ac- 
crete most of its mass. Very soon a disk and a bipolar outflow 
develop. The spectral energy distribution slowly shifts towards 
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shorter wavelengths: A class I protostar is already seen in the 
near-infrared even though its luminosity still is mostly due to 
accretion (at a canonical age of about 10 5 years). In the subse- 
quent classical and weak-line T Tauri stages (class II and III), 
most of the material has already been accreted leaving less and 
less circumstellar material. Flaring in YSOs is a ubiquitous phe- 
nomenon at diverse wavelengths, however it has not previously 
been studied simultaneously across the electromagnetic spec- 
trum in protostars and rarely in more evolved YSOs, even though 
multi-wavelength correlations would contain additional infor- 
mation on underlying physics. 

One would expect some correlation of thermal hot-plasma 
X-ray emission and non-thermal centimetric radio emission in- 
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terpreted as gyrosynchrotron radiation due to magnetic fields 
because both types of radiation are partially produced in the 
same regions. For a discussion of radio a nd X-ray obs erva- 
tions of (proto-)stellar coronae, see Gudei (l2002al l2004h. an d 
for a discussion of multi-wavelength aspects, see lGiide l (2002b). 
Looking at the sun as the star which is best studied for multi- 
wavelength variability, we see much uncorrelated variability but 
also a particular kind of radio-X-ray correlation in the form 
of the Neupert effect which was observed towards a star other 
than the sun (the dMe flare star UV Cet) for the first time by 
iGiidel et"aTI (fl996). Radio and X-ray emission both trace the in- 
nermost surroundings of protostars at less than a few stellar radii 
while near-infrared emission is mainly due to warm dust. 

A correlation of optical/near-infrared flares interpreted as ac- 
cretion events and X-ray variability would permit speculations 
on the connection of accretion and X-ray emission. Such a cor- 
relation can be expected from the observation of different X- 
ray luminosities for accreting and non-accreting T Tauri stars as 
well as a few cases where soft X-ray emission was interpreted 
as due to accretion (see Section jjl.ll i. Another cause for corre- 
lated variability is large cool star spots due to magnetic activ- 
ity. They lead to rotational modulation of photospheric light and 
may directly or indirectly influence X-ray emission. In the near- 
infrared, variations in the circumstellar disk emissio n are more 
impor tant than variability in the stellar photosphere ( Eiro a et all 
2002}. In case of X-ray emission due to accretion shocks (see 
next subsection), one would expect to see a correlation at near- 
infrared or even optical wavelengths, depending on the accretion 
rate. For the earliest evolutionary stages of stars, a direct com- 
parison of optical and near-infrared variability is not possible 
since they can only be observed in the near-infrared if at all. 
Variability at near-infrared and optical wavelengths can also be 
due to pulsation and eclipsing multiple systems. 

1.1. Radio and X-ray emission from YSOs 

YSOs are known to emit thermal as well as non-thermal cen- 
timetric radio emission. While the thermal emission appears 
to be mostly due to shock-induced ionisation at the base of 
outflows, the non-thermal emission mainly is gyrosynchrotron 
emission from mildly relativistic electrons gyrating in ma gnetic 
fields. For a review of radio properties of YSOs, see I Andrei 
( 1996 ); a review of associated hig h-energy processes was given 
bv lFeigelson & Montmerldd 1999ft . 

The main possibilities for the origin of the strong X-ray ra- 
diation from YSOs are (1) emission from a magnetically con- 
fined plasma (either in stellar coronae similar to that of the Sun 
and other late-type stars, or perhaps in magnetic structures con- 
necting the star and its circumstellar disk) with magnetic re- 
connection causing flares and (2) emission from hot accretion 
shocks at the stellar surface. Recent sensitive X-ray observa- 
tions of large samples of T Tauri stars (e. g. the Chandra Orion 
Ultradeep Project; see lGetman et a l. 2005) have now clearly es- 
tablished that the bulk of the X-ray emission from T Tauri stars 
must arise from magnetically confined plasma, most likely in 
rather compact coronal loops. Accretion shocks seem to con- 
tribute some fraction of the soft X-ray em ission of some TTS 
(e.g. TW Hya, XZ Tau and BP Tau, see iKastneret al.l 120021: 
Sch mitt et al.ll2.005h . but clearly cannot account for the high X- 
ray l uminosities and plasm a temperatures seen in the majority of 
TTS (Preibisch et al. 2005). Nevertheless, there is a relation be- 
tween X-ray activity and disk accretion in the sense that actively 
accreting TTS are weaker X-ray emitters than non-accreting 
TTS. This effect may be even more important for class I proto- 



stars. Possible explanations for this su ppression of X-ray emis- 
si on in accreting TTS a re discussed in Preib isch et al.l d2005l) and 
in lJardine et al.l (120061) . 

X-ray emission from class I protostars i s also well estab - 
lished observationally dlmanishi et al.l 120011; iPreibischl |2004), 
and an X-ray fl are towards a candidate class protostar in CrA 
was reported by Hamag uchi et al.l (120051) . The very high plasma 
temperatures derived from the protostellar X-ray spectra again 
clearly suggest an origin related to some kind of magnetic activ- 
ity, whereas accretion processes cannot reproduce the observed 
X-ray properties. 

Until now, there have been few simultaneous multi- 
wavelength observations of YSOs, mai nly targeting T Tauri stars 
and not class or class I protostars: Feigelson et al. ( 1994) ob- 
served a magnetically active T Tauri star (V773 Tau) simul- 
taneo usly at X-ray, ultraviol et, o ptical, and radio wa velengths, 
while iGullbring et al.l d 19971) and lStelzer etal.1 d2003l) study si- 
multaneous optical and X-ray variability of the classical T Tauri 
star B P Tau and the weak-li ne T Tauri star V410 Tau, respec- 
tively. Guenth er et al.l (HOOO) targeted pre-main sequence stars 
in the Taurus-Auriga region, specifically V773 Tau, in quasi- 
simultaneous X-ray, optical, and radio observations. In all cases, 
no correlated variability was fo und, based on relatively short 
observations. iBower et al.l d2003l) report the simultaneous radio 
and X-ray observation of an enormous radio flare towards a 
supposed highly obscured weak-line T Tauri star (GMR-A) in 
Orion. The X-ray flux increased by a factor of 8 two days be- 
fore the radio flare detection. The p Oph cloud complex and 
its pre-main sequence stars w ere selected for si multaneous ra- 
dio and X-ray observations bv lGagne et al.l d2004l) . covering the 
largest number of sources yet. However, none of th e class I pro- 
tostars were detected in both wavelength regimes. iBerger et al.l 
(2005) attempted to constrain the magnetic properties of a candi- 
date brown dwarf by simultaneous radio, X-ray, and Ho- obser- 
vations. Because of the low upper limit to the X-ray flux, they 
found a large radio -to-X-ray luminosity ratio. Most recently, 
Stassun et al. (2006) present a comprehensive optical variability 
study of Orion pre-main sequence stars, carried out simultane- 
ously with the Chandra Orion Ultra-Deep Project (COUP). They 
find indications for time-correlated optical and X-ray variability 
in only 5% of their sources and conclude that there is little evi- 
dence for a causal link between variability at optical and X-ray 
wavelengths. 

In a preparatory study to this multi- wave length project, we 
singled out the Corona Australis star-forming region, particu- 
larly the Coronet cluster, as the ideal target for this kind of ob- 
servation: It is a compact cluster of protostars and other YSOs, 
most of them detecte d at radio, infrared and X-ray wavelengths 
dForbrich et alj|2006t for a review of the entire Corona Australis 
star-forming region, see lNeuhauser & Fo rbrich 20071). 



2. Observations and data reduction 

We present here the results of a multi-wavelength observing 
campaign targeting YSOs in Corona Australis, based on five 
epochs of Chandra observations. For the first time, protostars 
(class 0/1 YSOs) were observed simultaneously at radio, X-ray, 
and near-infrared wavelengths. The field of view of the observa- 
tions is shown in Fig. Q] their scheduling is shown in Fig. [2] as 
well as listed in Table Q] 
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Fig. 1. Plot of the Coronet region cov- 
ered in this multi-wavelength dataset. The 
background image shows the combined X- 
ray data, the big circle indicates the half- 
power beam width of the VLA primary 
beam. The outer two boxes mark the cover- 
age of the SMARTS- UB VRI and SIRIUS- 
JHKs data (thick line). The two small 
boxes show the fields of view of the Soft 
observations. The central Coronet cluster 
is shown in detail in Fig. [3] Three 2MASS 
sources discussed in the text are marked 
as small circles: a) 2MASS J19014041- 
3651422 b) 2MASS J19013385-3657448 
c) 2MASS J19012872-3659317 



Due to the far southern declination of the Coronet clus- 
ter (nearly —37°), the region is accessible for the NRACQ Very 
Large Array only during a few hours per day. Thus, the core ex- 
periment was designed around five 4h duration Chandra obser- 
vations on successive days (with a one-day gap due to perigee 
passage when the satellite is inoperative). The instrument used 
was the Advanced CCD Imaging Spectrometer (ACIS). The 
Chandra observations were to be complemented by five epochs 
of simultaneous VLA observations around the Coronet Cluster's 
meridian passage. To ease scheduling, blocks of 2h VLA time 
were to be scheduled within the 4h Chandra observations. A 
few simultaneous near-infrared observations were collected with 
ESO's 3.5m New Technology Telescope (NTT) in Chile, using 
the SofI imager in the JHKs bands subsequently and in two 
fields east and west of R CrA. SofI is based on a 1024x1024 
Hawaii HgCdTe CCD detector, providing a field of view of 
4'94 x 4'.94- with a pixel scale of 0'.'288 (using the Large Field 
Objective). The exposure time was set to the minimal possible 
value at 1.182 seconds. 

Additional near-simultaneous observations were carried out 
with SIRIUS (Simultaneous Infrared Imager for Unbiased 
Survey) at the 1 .4m telescope of the InfraRed Survey Facility^ 
(IRSF) in South Africa (fo r a description of this simultaneous 
three-band NIR camera, see lNagavama et al.l2003h . The SIRIUS 
camera has a field of view of 7'7 xTJ and a pixel scale of 0'.'45. 
Four different fields were observed, one centered on the Coronet 
cluster, and three more north, west, and south of the central field. 
The SIRIUS images typically consist of ten dithered exposures 
of 30 sec duration. 



1 The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 

2 IRSF is operated by Nagoya University and SAAO (South African 
Astronomical Observatory) at Sutherland, South Africa. 



Table 1. Simultaneous observations of the Coronet with 
Chandra and the VLA (given in UT) 



Day(MJD') Chandra VLA 



Aug 


08, 2005 (0) 


02:37:52 


-07:33:11 


04:04:20 


- 06:03:45 


Aug 


09, 2005 (1) 


02:38:52 


-07:22:10 


04:00:50 


- 06:00:00 


Aug 


10, 2005 (2) 


01:58:25 


-06:53:28 






Aug 


12, 2005 (4) 


03:12:57 


-07:54:57 


04:18:30 


-06:18:20 


Aug 


13, 2005 (5) 


01:51:20 


-06:34:32 


03:44:40 


-05:44:10 



1 MJD-53590 



Lastly, simultaneous optical imaging (bands UBVRI succes- 
sively) of the region was performed with the CTIO 0.9m tele- 
scope, operated by the SMARTS consortium, at the Cerro Tololo 
Inter-American Observatory in Chile. The detector used was a 
2048x2046 CCD camera providing a field of view of 13:5 x 13:5 
with a pixel scale of 0V396. Again, two fields were observed, 
east and west of R CrA, but both overlapping to cover R CrA. 
The last three epochs including the last multi-wavelength epoch 
were increasingly compromised by bad weather. In all filters and 
both fields, long and short exposures were taken, e.g. for the I- 
band typically durations of 20 sec and 2 sec. 

2.1. X-ray observations 

The Chandra Interactive Analysis of Observations (CIAO) 3.3 
software package together with CALDB 3.2.1 was used for the 
analysis of the Chandra-ACIS data. Apertures of 5" in diameter 
were used for light curves and spectra alike. The spectra were 
compiled with the task psextract which includes the creation of 
response matrices and ancillary response files. Spectral bin sizes 
were set to ensure a minimum of 15 events per bin. 
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Fig. 2. Graphical depiction of the observing times at the different 
observatories indicated on the lefthand side. 



epoch, were used to establish an absolute flux density scale. The 
nearby phase calibrator 1924-292 was observed every seven 
minutes. We estimate that our absolute flux density scale thus 
derived is accurate to ~ 15 %. 

Due to the low elevation (< 19°), there are problems with 
phase instabilities caused by the atmosphere. After initial flag- 
ging of obviously faulty uv data, data were phase-only self- 
calibrated before imaging in Stokes-/ and Stokes- V could be car- 
ried out for different time intervals. Source positions and fluxes 
were determined by fitting Gaussians to each detection, using 
the MPS tasks SAD and JMFIT. Splitting the single observing 
epochs into intervals of approximately one hour, the root mean 
square (rms) noise level in the Stokes / and V maps is « 40 yuJy; 
for the entire uv data this value reduces to 17 /iJy. The synthe- 
sized beam size for the combined four epochs is 6 '.'7 x 2'.'0 at 
position angle 1°9. 
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Fig. 3. Combined Chandra data (Obs ID 5402-5406, smoothed). 
The circle denotes the approximate VLA half-power primary 
beam. 



2.2. Radio observations 

The VLA was observing in its D configuration (program ID S6 
255) at v = 8.44 GHz with two intermediate frequency (IF) 
pairs, offset by the IF bandwidth of 50 MHz. One IF of each 
pair detected right circular polarization (RCP), the other one left 
circular polarization (LCP). The phase center was (a, 5)j2ooo = 
19 h 01 m 48 s , -36°58'00". The primary beam of the VLA anten- 
nas limits the field of view to approximately five arcminutes 
FWHM around the phase center position (Fig. |3). Compared 
to the other observations reported here, the radio observations 
cover the smallest field, centered on the Coronet cluster. Data 
were analyzed using the NRAO Astronomical Image Processing 
System (AIPS). Observations of 3C286, once per observing 



2.3. Near-infrared and optical observations, photometry 

The SIRIUS-y/Z^Cs data were processed using a reduction 
pipeline developed as an external package of the IRAF0 soft- 
ware package. Processing includes dark frame subtraction, flat- 
fielding, sky subtraction and finally combination of the dithered 
single images. 

The Sofl-JHKs images, consisting of pairs of offset images, 
were processed using Eclipse 5.0.0. Dark frame subtraction, flat- 
fielding, and sky subtraction were performed. 

The optical images in U BVRI bands were dark-subtracted 
and flat-fielded using IRAF. Slightly different properties of the 
four CCD chips in the array had to be accounted for using the 
quadproc procedures. The observations were partly affected by 
bad weather in the last three epochs. 

By correlating the near-infrared and optical data with 
2MASS positions, an overall astrometric accuracy of < O'.'l was 
achieved. 

Differential photometry was carried out for optical and near- 
infrared co unterparts of X-ray sources using an algorithm de- 
scribed by iBroeg et all d2005l) . This procedure uses aperture 
photometry determined by the IRAF daophot package. For the 
aperture radii we chose values of around three times the radius 
of the point-spread function while the background was deter- 
mined from surrounding annuli. The respective comparison stars 
are weighted down according to their variability as represented 
by their standard deviations, resulting in an optimum artificial 
comparison star. As an estimate of the errors in the resulting 
light curves, a linear model of the instrumental errors is self- 
consistently adjusted so that it reproduces the standard devia- 
tions of the least variable comparison stars before it is applied 
to the science source. For the sources analyzed here, measure- 
ments with resulting lcr errors of more than 0.05 mag were dis- 
carded. Due to the different selections of comparison stars, the 
zero levels of the differential magnitude scales generally are not 
the same. Thus, the light curves are plotted with the respective 
median magnitudes subtracted. Sources observed outside the re- 
spective linear detector ranges were not analyzed photometri- 
cally. 



3 IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc. (AURA) under cooperative agreement 
with the National Science Foundation. 
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2.3.1. S CrA 

We undertook a separate analysis of the S CrA optical photom- 
etry in the context of our long term monitoring, which consists 
of over 200 observations between 2003 April and 2006 June, 
mostly obtained with the CTIO/S MARTS 1.3m ANDICAM 
dual-c hannel photometer. T his detector is described in, for ex- 
ample, IWalter et al.1 (120041) . Out of this long-term monitoring, 
we include seven ANDICAM data points in our analysis here, 
obtained in the nights of 2005 August 14-17. We do not make 
use of the short exposures because the data quality is poorer, and 
those images do not extend the temporal coverage. 

S CrA is spatially resolvable under the typical observing con- 
ditions at CTIO. We measure the brightness of the individual 
components by fitting two Gaussians to the image. The separa- 
tion and position angle are well-known, and are not free to vary. 
We constrain the two images to have the same point-spread func- 
tion (PSF). The relative brightnesses of the two stars are then the 
relative amplitudes of the two Gaussians. 

The seeing (measured as the FWHM of the Gaussian PSF) 
in the I band varied from 0'.'7 to 1'.'4, with a median value of 
0'.'9. The seeing was stable for the first 4 nights, with a variance 
of 0'.'06, but went up in the remainder of the run, with the vari- 
ance in the mean seeing ballooning to 0'.'2. There are no obvious 
seeing-dependent effects on the photometry. 

We calibrate the brightness of the system in a three-step pro- 
cess. First, we determine the long-term brightness variations of 
the S CrA binary through differential photometry with respect to 
9 other stars in the sparse field. We use an 8"aperture to enclose 
all the light from the S CrA binary. We assume that the weighted 
mean brightness of these nine stars does not vary significantly 
with time. The comparison brightness is heavily skewed to the 
flux of VSST 14, a pre-main sequence star, but not a known vari- 
able star, which is by far the brightest of the comparison stars in 
the field. To check whether our assumption of constancy is valid, 
we also compared the differential brightness against four other 
fainter stars singly. All yield the same long-term trends, suggest- 
ing that VSST 14 is not nearly as variable as S CrA and the 
assumption of constancy is warranted. 

We then use the mean long-term magnitude of S CrA 
(7=9.81) from the Van Vleck Observatory T Tauri databaseQ 
to establish the mean I-band magnitude. We assume that the 
mean differential magnitude in our photometric monitoring pro- 
gram corresponds to the mean long-term magnitude. We deter- 
mine this offset and apply it to all the differential measurements. 
Finally, we divide the light from S CrA between the two compo- 
nents based on the Gaussian fits. 



3. Single-wavelength variability 

After a discussion of the results by wavelength, we discuss 
multi-wavelength aspects for individual sources in section [4] 

3.1. Radio data 

The radio data show the same sou rces as in previous observa- 
tions (Fig. 2 in iForbrich et al.ll2006l) . although the D array con- 
figuration of the VLA used this time results in lower angular res- 
olution which is problematic only for IRS 7 and its immediate 
vicinity. The combined 2005 data is less sensitive than the com- 
bined 1998 data. No clear variability beyond an estimated uncer- 
tainty of ~ 15% in the absolute flux density scale is seen within 

4 http://www.astro.wesleyan.edu/~bill/research/ttauri.html 



Table 2. Radio flux densities of the Coronet sources 



Source ID 


2005 flux density 1 


1998 flux density- 




[mjy] 


[mJy] 


IRS 2 


0.49 ± 0.02 


0.29-0.36 


Source 5 


1.14 ± 0.02 


0.40-1.21 


IRS 5N 


0.06 ± 0.02 


0.10 3 


IRS 5 


0.68 ± 0.02 


0.75-3.28 


IRS 6 


0.07 ± 0.02 


0.15 


IRS 1 


0.42 + 0.02 


0.36-0.59 


R CrA 


0.30 ± 0.02 


0.16-0.26 


IRS 7W 4 


5.33 ± 0.02 


4.11-5.29 


IRS 7E 


1.77 ± 0.02 


1.38-2.10 



1 quoted errors are from Gaussian fitting. Our flux scale has an absolute 
uncertainty of ~ 15 % 

2 minima and max ima (if source detected in single epochs), see 
IForbrich et ail d2006h 

3 note that this source is not called IRS 5N in lForbrich et al.1 d2006h but 
is no. 5 in the tables 

4 blended with Source 9 



Table 3. X-ray count rates for the Coronet sources, from the 
five Chandra observations, 0.5-10 keV, corrected for effective 
exposure 



Source ID 


SI 


S2 


S3 


S4 


S5 




1/ksec 


1/ksec 


1/ksec 


1/ksec 


1/ksec 


IRS 2 


20.4 


27.2 


47.2 


31.3 


36.6 


IRS 5N 


0.1 


1.0 


0.8 


0.8 


0.0 


IRS 5 


16.1 


27.0 


49.0 


17.1 


11.7 


IRS 6 


1.3 


2.4 


2.2 


2.7 


2.0 


IRS 1 


25.3 


14.5 


14.5 


18.3 


26.7 


IRS 9 


0.5 


2.5 


0.7 


2.3 


0.9 


RCrA 


4.7 


6.4 


5.3 


23.7 


15.2 


IRS 7W 


1.0 


6.1 


5.1 


1.6 


1.9 


IRS 7E 


0.7 


0.7 


0.4 


1.3 


1.5 



these four epochs of observing which is why we only quote flux 
levels determined from the combined dataset in Table [2] while 
showing exemplary radio flux density curves for the three class I 
protostars as well as R CrA in Fig.|U Also on timescales shorter 
than the duration of each epoch (2h), no variability was found. 
For comparison, the minimum and maximum flux densities from 
the 1998 monitoring are shown as well in Table [2] The most 
striking differences are IRS 5 - the only known class I protostar 
with circularly polarized non-thermal centimetric radio emission 
- showing a lower flux than before throughout this observation 
and R CrA with a higher radio flux. The flux of IRS 5 is in fact so 
low that we cannot detect Stokes- V emission which presumably 
is below the detection threshold, which would indeed be the case 
if the polarization degree was comparable to the values measured 
in 1998. The only X-ray flare of IRS 5 during this observing 
run took place during VLA maintenance time so that it could 
not be monitored simultaneously at radio wavelengths. R CrA 
was marginally detected as a weak radio source before, now it 
is brighter and clearly detected. IRS 2 appears brighter than in 
the 19 98 data, similar to levels reported before (Feig elson et al.l 
119981) . When comparing the radio variability, it is important to 
note the different timescales probed in our 1998 and 2005 moni- 
toring campaigns. While in 1998, observations span nearly four 
months, in 2005 only five days are covered. Our results suggest 
generally lower levels of variability on these shorter timescales. 
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Fig. 4. Exemplary VLA radio flux density curves for the three 
class I protostars IRS 1, 2, and 5, as well as for R CrA. Flux 
densities were derived by fitting Gaussians with the AIPS task 
SAD. The error bars mainly take into account actual noise (lcr) 
rather than the quality of the fit. 



I 




RS 5 A/E 



Fig. 5. The double source IRS 5 (marginally resolved here) 
and its neighbouring source IRS 5N seen in Ks (Sofl/NTT) on 
August 9, 2005. 



3.2. X-ray data 

Focusing on multi-wavelength variability here, a detailed X- 
ray census of the Corona Australis region is beyond the 
scope of this study and will appear in a separate paper 
(Forb rich & Pre ibisch 2007). Compared to the previous data pre- 
sented by [Forbrich et al. (2006), several new X-ray sources have 
been detected in th e region of the Coronet cluster (Fig. 2 in 
Forb rich et al.l l2006 vs. Fig. [3]). With the exception of the source 
north-east of IRS 5 (IRS 5N, detected only in S2), none of these 
sources is detected simultaneously at radio, near-infrared or op- 
tic al wavelengths. IRS5N was detected as a weak radio source 
by iForbrich et al.l (|2006) and is also seen in the near infrared 
(Fig.0. 

Averaged count rates per observing epoch for the main 
Coronet sources are summarized in Table [3] T he count rates 
are si milar to previous Chandra count rates (e.g. Forbrich et al. 
2006). The most variable sources by factors of about five are 
IRS 5, R CrA, and IRS 7W. IRS 5 appears to be in a decay- 
ing flare state in epoch S3. R CrA shows increased levels of 
X-ray emission in the last two epochs S4 and S5. No strong 
flare towards IRS 7E occured during this observing run, but the 
so urce is detected with a lo w and slightly variable flux level (but 
cf. lHamaguchi et al.|[2005l) . The class I source IRS 9 shows an 
X-ray flare in the second observing epoch S2. Both IRS 7E and 
IRS 9 are too weak for the analysis of meaningful light curves. 
Thus we only give t he averaged count r ates. Confirming the de- 
tection reported by Ski nner et alj d2004l) . we also see T CrA in 
X-rays, south-east of IRS 7. We note that the s pectral proper- 
ties of the sources are similar to those reported by Forbrich et al. 
(2001). 



There are a few interesting sources outside the VLA primary 
beam which merit attention in this context because they are vari- 
able and were observed in at least two wavelength regimes si- 



multaneously: S CrA shows variability in the first half of the 
observations and 2MASS J19012872-3659317 shows a flare in 
the last Chandra observing epoch S5. The X-ray light curves are 
plotted in Figs. [7] and [8] 

3.3. Near-infrared and optical data 

In a first step, near-infrared and optical data were searched for 
counterparts of X-ray sources. We then singled out those sources 
that were variable by more than 0. 1 mag in the respective band. 
The discussion focuses on these variable sources. We note, how- 
ever, that while a full discussion of all sources is beyond the 
scope of this paper, not all X-ray sources in this region have NIR 
or optical counterparts and vice versa. In general, only slow vari- 
ability is seen at these wavelengths, with the notable exception 
of TY CrA which showed an eclipse minimum. 

Again, many sources are not covered by the VLA primary 
beam, but there is a good near-infrared/optical to X-ray corre- 
spondence in source coverage. Fifteen near-infrared counterparts 
of X-ray sources spread over the four different pointings of the 
SIRIUS imager were analyzed for variability in the J band, spar- 
ing seven sources which were too bright in the J band for analy- 
sis in the / band (see below, a few of them were in fact also too 
bright in the / band). Two class I protostars, namely IRS 5 and 
IRS 9, as well as the embedded source north-east of IRS 5, are 
undetected in J band and thus had to be analyzed in Kg band. 

Among the selected SIRIUS /-band sources, four show vari- 
ability of > 0. 1 mag (see Fig|7]i, namely the two class I protostars 
IRS 1 and IRS 2 as well as 2MASS J19013385-3657448, and 
2MASS J19014041-3651422. It is remarkable that this selection 
includes the two identified class I protostars which are detected 
at all in the J band. In the K s band, IRS 9, IRS 5N and IRS 5 
(the latter being shown in Fig. [6j vary by less than 0.1 mag, but 
IRS 9 and IRS 5N brightened steadily during the observations. 
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Fig. 6. Top panel: Chandra X-ray light curves 
of IRS 5 for the five observations S1-S5, bin 
size is 1 ksec. Center and bottom panel: Ks- 
band differential photometry of IRS 5 taken 
with SIRIUS and Sofl, respectively. If any- 
thing, the K s -band flux appears to slightly de- 
crease around the enhanced X-ray activity of 
IRS 5. Error bars are lcr. 



Most of the saturated /-band counterparts of X-ray sources 
except for HD 176386 (mv=7.30) were also detected in / band 
and available for photometry. Several sources show variability 
of > 0.1 mag, namely the Herbig Ae stars R CrA, T CrA, and 
TY CrA, as well as the T Tauri star S CrA and the embed- 
ded sources ISO-CrA 1 37 and HBC 680, an M3 dwarf WTTS 
dNeuhauser et al.l l2000b . The light curves are shown in Fig. [8] 
and in Fig. [HI 



4. Multi-wavelength variability 

4.1. The class I protostars IRS 2,5, 1, and 9 

The purpose of these observations was to, for the first time, ob- 
tain simultaneous multi-wavelength data on class I protostars. As 
stated above, no significant radio variability was observed dur- 
ing the four VLA observations. IRS 2 and IRS 1 show moderate 
X-ray variability (see Fig. |7), together with uncorrelated near- 
infrared variability, while IRS 5 is again the most variable class I 
protostar in X-rays. The enhanced activity in Chandra observa- 
tion S3 was not covered by simultaneous radio observations. The 
Ks -band flux of IRS 5 is slightly lower before and during the 
S3 observation (Fig. |6). The radio luminosity of IRS 5 is much 
lower during the observations reported here than before, with its 
Stokes- V emission below the detection limit. IRS 5N, the weak 
X-ray source north-east of IRS 5 (Fig. ID is detected in X-rays 
only in epoch S2. The source was simultaneously observed in the 
near-infrared with Sofl/NTT, however no significant variability 
could be found compared to the Sofl/NTT exposures taken later 
on except for a slight steady brightening by < 0. 1 mag through- 
out the SIRIUS ^-band observations. While the class I source 
IRS 9 shows an X-ray flare in the S2 Chandra observation, it re- 
mains undetected in the VLA observations as before. No signifi- 
cant near-infrared variability besides a slight steady brightening 
by < 0.1 mag is seen in the SIRIUS Ks -band data (the X-ray 
flare was not covered simultaneously at NIR wavelengths, both 
light curves not shown). 

4.2. The Herbig Ae/Be stars R CrA, T CrA, and TY CrA 

The three Herbig stars R CrA (type Ae), T CrA (type FOe), 
and TY CrA (type R8e) are among the most variable /-band 
sources in this dataset. R CrA which appeared t o be in differ- 
ent X- ray activity stages in the data presented by iForbrich et alJ 
(2006) is flaring continuously and is particularly active in the 
last two Chandra (S4-5) observations (Fig. [8]). The X-ray emis- 
sion of R CrA is again associated with an extremely high plasma 
temperature. Simultaneously, R CrA continuously brightened 
by approximately 0.6 mag (/ band) throughout the observation 



campaign. The optical data, however, are compromised by bad 
weather directly after the S4 Chandra observation. Still it seems 
that the brightening stops after the onset of enhanced X-ray 
activity, if only to resume later. The same brightening is ob- 
served in the U band (not shown). Given the earlie r discussion on 
the or igin of the X-ray emission towards R CrA (Forbric h et al.l 
2006), it is quite likely that two different objects are in fact ob- 
served here in optical and X-ray emission. 

T CrA is quite variable at optical wavelengths (Fig.lTTI). how- 
ever it is only barely detected as an X-ray source. While our ob- 
serving run spans several eclipse minima of the multiple Herbig 
Be system TY CrA, no minimum has been observed in different 
bands simultaneously (see optical data in Fig. [8). 

4.3. The classical T Tauri stars S CrA and IRS 6 

S CrA is a 1 '.'3 visual pair, varying by up to two magnitudes and 
showing variable inverse P Cygni profiles in the upper B almer 
lines, indicative of mass accretion. iPrato & Simonld 19971) found 
that both components can be identified as TTSs. The source re- 
mains unresolved in our X-ray data due to its off-axis location. 
Our X-ray light curve shows a quite slowly varying count rate 
with an apparently uncorrelated /-band light curve (see Fig . [HJ . 
The shape of the optical light curve is quite complex though not 
unusual. 

The /-band light curves of the NW and SE components are 
shown in comparison to the X-ray data Fig.[8]while Fig.|9]con- 
tains the full light curve discussed here together with the mag- 
nitude difference of the two components. Both components of 
the binary brightened throughout the time of the X-ray obser- 
vations, before drastically fading. At peak the system is about 
0.3 mag brighter than the historical mean, and within 0. 1 mag of 
the maximum brightness we have observed. There is a 0.5 mag 
brightening in the NW component with a duration < 2 days, 
on the second day. The NW component is consistently brighter 
than the more volatile SE component. Historically (that is, in 
our database), the NW component is brighter than the SE com- 
ponent at / 89+2% of the time, with a median brightness dif- 
ference of 0.45 mag. We observed a median brightness differ- 
ence of 0.48 mag. The rate of change of the stars is also well 
within the norms. The only oddity is the slow apparently coher- 
ent variations in the two stars. We would be concerned if the 
stars varied exactly in phase, but Fig. |9] shows that throughout 
this time the SE component is brightening relative to the NW 
component. Still, the concordance of the component light curves 
is unusual. The light travel time between the components, at a 
projected separation of 170 AU, is about one day. This makes it 
difficult to explain the rapid fading in terms of, say, a dust event 
in the system. 
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Fig. 7. Chandra X-ray light curves for 
YSOs with near-simultaneous IRSF J- 
band data in the respective lower panels 
(with lcr error bars). 



Even in this more detailed photometric analysis of the two in the NW component, there does not seem to be a direct con- 
components of S CrA, we find no apparent correlation between nection (see Fig.fTOli. 

the X-ray light curve and the optical light curve. While the , , . ,. 

, ■ , \r ■ , , , , „ , a The T Taun star IRS 6 is detected in centimetric radio emis- 

bnghtest X-ray points correspond closely to the 0.5 mag flare . ... , , , e , ■ , . , i onn A t i i 

b J r r J b sion, at a similar level as before Forbrich et al. 2006). Low-level 



J. Forbrich et al.: Simultaneous multi-wavelength observations of YSOs in the Coronet cluster 
a) R CrA 



9 



■ 0.04 r 
0.03 - 

o.oa - 



§ 0.02 - I 



- 
i.5 - 



0.0S 7 

' 0.015 ~ 7 

0.01 - 

0.005 ~ 7 

L 

-0.05 1 

- 

0.05 1 + 

0.1 - 



MJD 53590 

b) S CrA 
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c) TY CrA 



MJD 53590 

d) ISO-CrA 137 



Fig. 8. Chandra X-ray light curves for 
YSOs with simultaneous CTIO /-band 
data in the respective lower panels (with 
lcr error bars). The /-band light curve for 
S CrA results from a special analysis in- 
volving Gaussian fitting to its two compo- 
nents NW and SE (see Fig. [9] and text). 
Figure continues on next page. 
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X-ray variability is accompanied by insignificant changes in the 4.4. Additional sources 
near-infrared emission during the observation run. 

Four more sources merit attention in this context because 
of peculiar X-ray or optical/NIR variability. The averaged 
X-ray count rates for these sources are listed in Table [4] 
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Fig.|8j Continued. 



2MASS J19013385-3657448 is observed at low and barely 
changing X-ray count rates, however shows a maximum J- 
band magnitude approximately in the middle of the observa- 
tions. 2MASS J19014041-3651422 drops continuously in J- 
band brightness, however shows fairly constant X-ray count 
rates. More variable X-ray emission with enhanced activity in 
epoch S3 is observed towards 2MASS J19012901-3701484 = 
ISO-CrA 137, a source showing a well-defined /-band emission 
maximum in-between Chandra observations S2 and S3. While 
2MASS J 190 12872-36593 17 shows the most pronounced X-ray 
variability of all these sources, the /-band light curve, while 
structured, shows only < 10% (0.1 mag) variations. Only one 
/-band measurement was taken during the enhanced X-ray ac- 
tivity, and this data point has a large error. The flare in epoch 
S5 is not covered in U band due to bad weather. The M dwarf 
HBC 680 is a rapidly variable source in the / band (Fig. fTTT i. 
but was not detected in X-rays, likely due to its proximity to the 
edge of the field-of-view. 
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5. Conclusions 

YSOs in and around the Coronet cluster in Corona Australis, 
including class I protostars, were observed simultaneously in 
multi-epoch X-ray, radio, near-infrared and optical observations. 
The main conclusions are as follows: 

- The sources in the Coronet cluster are variable in X-rays and 
near-infrared/optical emission, however no variability at ra- 
dio wavelengths could be detected even though the measured 
flux densities differ from previously measured values, indi- 
cating long-term variability. 

- No clear multi-wavelength correlations in variability were 
found. The near-infrared variability of the class I protostars 
observed here appears to be uncorrected to their X-ray vari- 
ability (with the exception of a slight - possibly insignificant 
- drop in Ks band flux during enhanced X-ray activity of 
IRS 5). This is what would be expected in case of coronal 
X-ray emission - as opposed to X-rays due to accretion - 
which is also indicated by the very high plasma temperatures 
involved. 

- Towards the Herbig Ae star R CrA, an apparent jump in X- 
ray emission is observed during continuous brightening in 
the optical regime. While the classical T Tauri star S CrA 
shows a complex optical light curve, its X-ray emission is 
much less variable. Its two components were found to nearly 
vary synchronously in the / band. 

- The radio behaviour of the polarized radio continuum source 
IRS 5, a class I protostar, could not be studied due to the 



Fig. 9. Upper panel: The total /-band magnitude of S CrA, deter- 
mined as described in the text. The dashed horizontal line is the 
historical mean in the Van Vleck Observatory data base, used as 
the zero point here. Lower panel: The magnitude difference be- 
tween the two components of S CrA. The dashed horizontal line 
is the median difference. 



low flux density of the source. The source shows, however, 
enhanced X-ray activity in one observation. 
- Among the four sources with X-ray counterparts and an near- 
infrared (J) variability exceeding 0.1 mag, there are two 
class I protostars, a disproportionate share in spite of the 
low numbers. In case of the / band variability, six sources 
with variability exceeding 0.1 mag include three Herbig 
Ae/Be stars, one CTTS, as well as two additional embedded 
sources. 
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Table 4. X-ray count rates for discussed sources outside the VLA primary beam, from the five Chandra observations, 0.5-10 keV, 
corrected for effective exposure 



Source ID 


SI 


S2 


S3 


S4 


S5 




1/ksec 


1/ksec 


1/ksec 


1/ksec 


1/ksec 


2MASS J19012872-3659317 


21.3 


28.5 
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26.5 


67.4 


2MASS J19012901-3701484 = ISO-CrA 137 
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2MASS J19013385-3657448 


3.3 
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Fig. 10. Detail of the X-ray and /-band light curves shown in 
Fig. [8] The optical flare seen in the NW component does not 
appear to directly correspond to the maximum in the X-ray light 
curve. 
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